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Q ■ We show that the NLS systems with multiplicative noise, nonlinear 

, ^, . damping and nonlocal dispersion exhibit a variety of interesting 

effects which may be useful for modelling the dynamical behavior 
K^ of one- and two-dimensional systems. 
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vn 
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O ' We consider the nonhnear Schrodinger (NLS) system with cubic non- 

linearity. This system models optical media, molecular thin films in the 
continuum limit, deep water waves and many other physical systems which 
^ i exhibit weak nonlinearity and strong dispersion. The effects we shall discuss 

are observed in our earlier works [^ ^ H and cited therein. 

'^ ■ The motivation for studying the 2-D NLS equation with thermal fluctua- 

tions is our intension to model efficient energy transfer in Scheibe-aggregates. 
O . Our starting point is a two-dimensional Davydov model with nonlinear cou- 

pling between the exciton and phonon system, and white noise in the phonon 
system. In Section 2 we derive a single equation for the exciton system with 



1 INTRODUCTION 



^ ■ multiplicative colored noise and a nonlinear damping term. In the continuum 

limit the collective coordinate approach indicates that an energy balance be- 
tween energy input (from the noise term) and dissipation can be established. 
Thus, this model may describe the state of thermal equilibrium in the molec- 

is 



ular aggregate. The coherent exciton moving on the aggregate [|T0|, |TT 



modelled by the ground state solution to the 2-D NLS equation, and the 
lifetime has been related to the collapse time of the ground state 0]. For 
sufficiently strong nonlinearity the thermal fluctuations will slow down the 
collapse. As a result this lifetime increases with the variance of the fluctua- 
tions, i.e. the temperature. 

In Section 3 the role of long-range dispersive interaction in the one- 
dimensional molecular system is investigated. Dispersive interactions of two 
types (the power and the exponential dependences of the interaction intensity 
on the distance) are studied. If the intaraction decreases with the distance 
slowly, there is an interval of bistability where two stable stationary states: 
narrow, pinned states and broad, mobile states exist at each value of the ex- 
citation energy. For cubic nonlinearity the bistability of the solitons occurs 
already for dipole-dipole dispersive interaction. We demonstrate a possibil- 
ity of the controlled switching between pinned and mobile states applying 
a spatially symmetric perturbation in the form of a parametric kick. The 
mechanism could be important for controlling energy storage and transport 
in molecular systems. 

2 NOISE AND DAMPING 

Following the derivation given in Ref. []1|, we start by assuming that the 
coupled exciton-phonon system can be described by the following pair of 
equations 

ihlpn + Yl Jnn'i^n' + XUni^n = , (l) 

n' 

Miln + MXiln + Mulun - x\'4'n? = Vn{t) ■ (2) 

Here ipn is the amplitude of the exciton wave function corresponding to site 
n and m„ represents the elastic degree of freedom at site n. Furthermore, 
— Jnn' is the dipole-dipole interaction energy, x is the exciton-phonon cou- 
pling constant, M is the molecular mass, A is the damping coefficient, ujq is 
the Einstein frequency of each oscillator, and rin{t) is an external force acting 
on the phonon system. To describe the interaction of the phonon system 
with a thermal reservoir at temperature T, ?7„(t) is assumed to be Gaussian 
white noise with zero mean and the autocorrelation function 

{nn{t)r^n'{t')) = 2MXkBT5{t - t')5nn' , (3) 



in accordance with the fluctuation-dissipation theorem ensuring thermal equi- 
librium. 

In order to derive a single equation for the dynamics of the exciton system, 
we start by writing the solution to Eq. (Q) in the integral form. Neglecting all 
exponentially decaying transient terms and making the additional assump- 
tion that tpn varies slowly in space and that only nearest-neighbor coupling J 
is of importance, we obtain in the continuum approximation for the con- 



tinuous exciton field il){x,y,t) = JV/Je '^*''*/^'?/^„(t) the equation of motion 

#t + V'V^ + l^l'V' - A^dV^n* + dV^ = , (4) 

where x and y are scaled on the distance i between nearest neighbors, a noise 
density a{x,y,t) is not white, but strongly colored |I|, and time was trans- 
formed into the dimensionless variable: Jt/h -^ t. The nonlinear damping 
parameter A = ■^. 

It can easily be shown that in spite of the presence of the nonlinear 
damping and multiplicative noise terms in Eq. (^), the norm, defined as 

N = J J \^{x,y,t)\^dxdy (5) 

will still be a conserved quantity. 

To investigate the influence on the collapse process of the damping and 
noise terms in Eq. (^, we will use the method of collective coordinates. 
To this end, we will make some simplifying assumptions. We will assume 
isotropy, which effectively reduces the problem to one space dimension with 
the radial coordinate r = \^x'^ + y'^. We also assume that the noise a can 
be approximated by radially isotropic Gaussian white noise. The validity of 
the approximation was discussed in Ref. |]I|. Finally, we assume that the 
collapse process can be described in terms of collective coordinates using the 
following self-similar trial function for the exciton wave function V'(r, t) 

^{r,t) = A{t)sech(-^]e^'^'^'>\ (6) 

This trial function, with three real time- dependent parameters A, B, and a 
determining the amplitude, width, resp. phase of the wave function, was used 



in Refs. |1l2| , g to investigate the case when A = in Eq. (§). The choice of 



this particular type of trial function can be motivated by regarding it as a 
generalization of the approximate ground state solution to the ordinary 2-D 
NLS found in Ref. |^. From the definition (^ of the norm, we immediately 
obtain the relation between amplitude and width A{t) ~ \fN/B{t). In 
analogy with the treatment for the undamped case in Ref. 0], we find that 
it is possible to arrive to the following differential equation for the width B 
of the exciton wave function 

53 B^ B^ ' ^ ' 

where the constants A and F are some functions of A^ and A. Note that 
A and F depend on the initial conditions via A^, and that while A can be 
either positive or negative, F is always positive. In the absence of noise and 
damping the collapse will occur if and only if A < 0. The white noise h{t) 
has the autocorrelation 

{h{t)h{t')) = 2D5(t-t') , (8) 

where D is the dimensionless noise variance. 

Our numerical calculations [^ show that for D < -Dp^it — 0.15, the effect 
of the noise is to delay the pseudo-collapse in terms of the ensemble average 



of the width, in analogy with the similar result obtained in Ref. [0] for the 
undamped case. For D > -Dgrit' ^^ observe a non-monotonic behavior of 
{B{t)). Initially, the average width will decrease in a similar way as when 
D < -Dctji;, but after some time {B{t)) will reach a minimum value and 
diverge as t ^ 00. This is due to the fact that for D > D^^^^, the noise is 
strong enough to destroy the pseudo-collapse and cause dispersion for some 
of the systems in the ensemble. As t — > 00 the dominating contribution 
to {B{t)) will come from the dispersing systems for which B ^ 00, and 
consequently {B{t)) will diverge for D > Df^j.^^. The minimum value of 
{B(t)) will increase towards -B(O) as D increases. 

3 NONLOCAL INTERACTIONS 

In the main part of the previous studies of the discrete NLS models the 
dispersive interaction was assumed to be short-ranged and a nearest-neighbor 
approximation was used. However, there exist physical situations that def- 
initely can not be described in the framework of this approximation. The 



DNA molecule contains charged groups, with long-range Coulomb interaction 
(1/r) between them. The excitation transfer in molecular crystals and the 
vibron energy transport in biopolymers are due to transition dipole-dipole in- 
teraction with 1/r^ dependence on the distance, r. The nonlocal (long-range) 
dispersive interaction in these systems provides the existence of additional 
length-scale: the radius of the dispersive interaction. We will show that it 
leads to the bifurcative properties of the system due to both the compe- 
tition between nonlinearity and dispersion, and the interplay of long-range 
interactions and lattice discreteness. 

In some approximation the equation of motion is the nonlocal discrete 
NLS equation of the form 

^-IT^n + Yl •Jn-m{i^m " ^n) + l^nl^^n = , (9) 

where the long-range dispersive coupling is taken to be either exponentially, 
Jn = Je~^l"l, or algebraically, Jn = J\n\~'^, decreasing with the distance 
n between lattice sites. The parameters /? and s are introduced to cover 
different physical situations from the nearest-neighbor approximation {f3 -^ 
oo, s — s> oo) to the quadrupole-quadrupole (s = 5) and dipole-dipole (s = 3) 
interactions. The equation (^ conserves the number of excitations A^ = 

n 

We are interested in stationary solutions of Eq. (^ of the form ipn{t) = 
0„exp(zAt) with a real shape function 0„ and a frequency A. This gives the 
governing equation for (pn 

A(j)n = E Jn-m{<Pm ' <Pn) + <PI ■ (10) 

Figure 1^ shows the dependence A^(A) obtained from direct numerical solution 
of Eq. (p!0| ) for algebraically decaying Jn-m- A monotonic function is obtained 
only for s > Scr- For 2 < s < Scr the dependence becomes nonmonotonic 
(of A/'-type) with a local maximum and a local minimum. These extrema 
coalesce at s = Scr — 3.03. For s < 2 the local maximum disappears. The 
dependence A^(A) obtained analytically using the variational approach is in 
a good qualitative agreement with the dependence obtained numerically (see 
0). Thus the main features of all discrete NLS models with dispersive in- 
teraction Jn-m decreasing faster than |n — m|~'*'='' coincide qualitatively with 




Figure 1: Number of excitations, A^, versus frequency, A, found numerically 
from Eq. ( [TI1| ) for s = oo (full), 4 (dotted), 3 (short-dashed), 2.5 (long- 
dashed), 2 (short-long-dashed), 1.9 (dashed-dotted). 



the features obtained in the nearest-neighbor approximation where only one 
stationary state exists for any number of excitations, A^. However in the 
case of long-range nonlocal NLS equation (|^), i.e. for 2 < s < Scr, there 
exist for each N in the interval [Ni{s) , Nu{s)] three stationary states with 
frequencies Ai(A^) < A2{N) < A^lN). In particular, this means that in the 
case of dipole-dipole interaction (s = 3) multiple solutions exist. It is note- 
worthy that similar results are also obtained for the dispersive interaction 
of the exponentially decaying form. In this case the bistability takes place 
for /3 < 1.67. According to the theorem which was proven in [^, the neces- 
sary and sufficient stability criterion for the stationary states is dN/dA > 0. 
Therefore, we can conclude that in the interval [iV;(s), iV„(s)] there are only 
two linearly stable stationary states: Ai(A^) and A3(A^). The intermediate 
state is unstable since dN/dA < at A = A2. 

The low frequency states are wide and continuum-like while the high fre- 



quency solutions represents intrinsically localized states with a width of a few 
lattice spacings. It can be shown that the existence of two so different soliton 
states for one value of the excitation number, N, is due to the presence of two 
different length scales in the system: the usual scale of the NLS model which 
is related to the competition between nonlinearity and dispersion (expressed 
in terms of the ratio N/J ) and the range of the dispersive interaction ^. 

Having established the existence of bistable stationary states in the non- 
local discrete NLS system, a natural question that arises concerns the role of 
these states in the full dynamics of the model. In particular, it is of interest 
to investigate the possibility of switching between the stable states under the 
influence of external perturbations, and to clear up what type of perturba- 
tions can be used to control the switching. Switching of this type is important 
for example in the description of nonlinear transport and storage of energy 
in biomolecules like the DNA, since a mobile continuum-like excitation can 
provide action at distance while the switching to a discrete, pinned state can 
facilitate the structural changes of the DNA 0. As it was shown recently 
in 0], switching will occur if the system is perturbed in a way so that an 
internal, spatially localized and symmetrical mode ('breathing mode') of the 
stationary state is excited above a threshold value. 
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